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Abstract. Bacterial flagella assume different helical shapes during the tumbling phase of a bacterium but 
also in response to varying environmental conditions. Force-extension measurements by Darnton and Berg 
explicitly demonstrate a transformation from the coiled to the normal helical state [N. C. Darnton and 
H. C. Berg, Biophys. J. 92, 2230 (2007)]. We here develop an elastic model for the flagellum based on 
Kirchhoff's theory of an elastic rod that describes such a polymorphic transformation and use resistive force 
theory to couple the flagellum to the aqueous environment. We present Brownian dynamics simulations 
that quantitatively reproduce the force-extension curves and study how the ratio F of torsional to bending 
rigidity and the extensional rate influence the response of the flagellum. An upper bound for F is given. 
Using clamped flagella, we show in an adiabatic approximation that the mean extension, where a local 
coiled-to-normal transition occurs first, depends on the logarithm of the extensional rate. 
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1 Introduction 

Many types of bacteria, such as Escherichia coli and Sal- 
. monella typhimurium, propel themselves forward by rotat- 
ing a bundle of elastic filaments with helical shape. Each of 
these flagella, as they are called, is driven by a rotary mo- 
tor. When the sense of rotation of one motor is reversed, 
the attached flagellum leaves the bundle and undergoes a 
sequence of different helical configurations characterized 
by their pitch, radius and helicity. During the flagellar 
polymorphism, the bacterium tumbles and then contin- 
ues swimming in a different direction, when the flagel- 
lum assumes its original form and returns into the bundle. 
Whenever the bacterium senses a positive food gradient, 
it prolongs the swimming phase and thereby increases the 
food uptake [1]. In recent years, bacteria and their flagella 
have also been used in studies relevant to microfluidics 
to transport colloids [2,3] and pump fluids [4] but also to 
create new liquid- crystalline phases of screw like objects 
[5]. An understanding of the bacterial tumbling motion 
and the applications of bacterial flagella in nanotechnol- 
ogy and microfluidics requires a sufficiently simple elastic 
model that includes the flagellar polymorphism. This ar- 
ticle aims to provide such a model. 

The bacterial flagellum consists of three parts. The ro- 
tary motor is embedded in the cell membrane and trans- 
mits its motor torque to the long helical filament with the 
help of the short and very flexible proximal hook. The 
filament of E.coli bacteria is up to 10 |ixn long and is 
about 0.02 (xm in diameter. It is relatively stiff but can 
switch between distinct polymorphic forms. The filament 



assumes these forms in response to external perturbations 
such as changes in pH value, salinity, and temperature of 
the surrounding fluid [6-8] , the addition of alcohols [9] or 
sugars [10], and by applying external forces or torques to 
the filament [11-14]. 

The helical filament is a cylinder formed by 1 1 protofil- 
aments each of which consists of a stack of protein mono- 
mers called flagellin. These monomers assume two con- 
formations (L and R) that mainly differ in length. Each 
protofilament only contains one type of monomer and there- 
fore L- and R-state protofilaments have different lengths. 1 
Mixing them in the flagellar filament produces bending 
which together with an intrinsic twist of the filament gives 
a helical configuration. In total two straight and 10 heli- 
cal polymorphic states are possible since protofilaments 
of one type cluster (see Fig. 1). Most of them were ob- 
served experimentally [15]. This is the picture Asakura 
[16] and later Calladine [17], with a more detailed geo- 
metric model, developed to explain the flagellar polymor- 
phism. The molecular structure for the flagellar filament 
was confirmed recently [18, 19] and several extensions of 
Calladine's model exist [15,20-23]. 

This article is mainly motivated by recent experiments 
of Darnton and Berg [14]. With the help of an optical 
tweezer set up they pulled the two ends of the flagellar fil- 
ament apart with a constant velocity and induced a tran- 



1 This is one important ingredient of the model developed 
by Asakura and Calledine. The different helical configurations 
following from this model were confirmed later for example in 
Ref. [18]. 
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Fig. 1. Twelve polymorphic states of the bacterial fiagcllum 
with curvature k — K; max sin(7rn_R/ll) and torsion r = tl + 
(tr — TL)nn/ll after Calladine [17], where na is the number 
of protofilaments in the R state. The quantities K max ~ 2.4/u.m, 
tl ~ — 5.2/u.m, and tr w 11.8/u.m are fit parameters. 



sition between two polymorphic configurations. They then 
reversed the velocity to compress the flagellum in order to 
return it to the initial configuration. Darnton and Berg 
recorded force-extension curves mainly for the transfor- 
mation from the coiled to the normal configuration. The 
transformation starts locally at one end of the flagellum 
and then proceeds in discrete steps along the flagellum. 
Signatures of the steps are clearly visible in the force- 
extension curves. 

Calculating the force-extension curves on the basis of 
coarse-grained molecular dynamics simulations is not pos- 
sible since simulation times are far below the experimen- 
tally relevant time scale of one second [24] . Modeling the 
polymorphism of a bacterial flagellum on a mesoscopic 
level is more appropriate. Goldstein et al. extended Kirch- 
hoff's classic theory of an elastic rod by introducing a 
double-well potential for the spontaneous torsion to de- 
scribe the transition between two helical states [25,26]. 
Wada and Netz also described the helical filament by Kirch- 
hoff's rod theory but attached a spin variable along the 
filament in order to distinguish locally between the two 
helical states [27] . They then performed hybrid Brownian- 
dynamics Monte-Carlo simulations to numerically calcu- 
late force-extension curves of bacterial flagella. 

In this article we present Brownian-dynamics simula- 
tions of the force-extension curves based on a model that 
is less time-consuming than the approach of Ref. [27] but 
that is completely equivalent as we demonstrate in an ap- 
pendix. Wc furthermore concentrate on different aspects 
of the force-extension curves, namely how they depend on 
the ratio of torsional and bending rigidities r and on the 
velocity or extensional rate with which the flagellum is 
pulled apart. We also give an upper bound for r which is 
partially in contrast to experimental results. The mean ex- 
tension, at which a coiled-to-normal transition first occurs 
locally, is a function of the extension rate. We demonstrate 
how this extension can be inferred from equilibrium prop- 
erties of a clamped helical filament. Our modeling is in the 
spirit of Refs. [25,26]. However, we show that a conven- 
tional double-well potential cannot reproduce the exper- 
imentally observed force extension curves. We therefore 
developed an alternative model, where we just "glue" the 




Fig. 2. The conformation of a slender elastic rod is described 
by the space curve r(s) of its center line and the material 
frame {ei,e2,ea}. The vector nn — d a e.z describes the local 
curvature of r(s). 



harmonic elastic energies of the two helical states together. 
We perform our simulations with realistic parameter val- 
ues and can directly compare our results to experiments. 

The article is organized as follows. In Section 2 we 
present our extension of Kirchhoff's clastic rod theory 
to include the polymorphism of helical filaments, explain 
how we perform the Brownian dynamics simulations, and 
present analytic results for a uniformly stretched conven- 
tional helical filament. Section 3 discusses the results for 
the simulated force-extension curves and addresses the 
clamped filament. We close with a summary and conclu- 
sions in Section 4. 



2 Continuum model for the bacterial 
flagellum 

2.1 Classic theory of an elastic rod 

The conformation of a slender rod with contour length 
L is described by the space curve of its center line r(s), 
parametrized by the arc length s, and a material frame 
of three orthogonal unit vectors {ei,e2,e 3 }, attached to 
each point on the center line. The vector points along 
the tangent of r(s) and ei, e 2 typically correspond to the 
principal axes of the inertia tensor of the cross section 
as indicated in Fig. 2. 2 Since the are unit vectors, the 
transport of the material frame along the center line is 
described by the generalized Frenet-Serret equations 



d s e.i — fl x 



(1) 



where d s means derivative with respect to s. So the con- 
formation of an clastic rod or filament is completely char- 
acterized by the angular strain vector fl = i?2, ^3). 
Alternatively, one uses the Frenet frame consisting of the 
tangent vector t = es, the curvature vector nn = d s t, and 
the binormal b = tx n. The Frenet frame transforms into 
the material frame when it is rotated about t = e$ by 



For a circular cross section the eigenvalues of the inertia 
tensor degenerate and we are free to choose the vectors e.\ and 
e.i- 
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the twist angle <fi, as indicated in Fig. 2. This relates the 
vector fl to the curvature k and torsion r of the filament: 

fli — ksiik/), i?2 = KCOS0, i?3 = t + d s 4>, (2) 

where the components fli are given with respect to the 
local material frame. 

Kirchhoff 's theory expands the elastic free energy T of 
a deformed elastic rod up to second order in the angular 
strain fl, 

F=l f cl (fl,fl )ds (3) 
Jo 

f cl (fl, fl ) = 4(^i) 2 + - ko) 2 + §(^ 3 - T ) 2 , 



where we also introduced a spontaneous curvature kq and 
torsion tq within S~1q = (0, Kq,tq), and A is the bending 
rigidity and C the torsional rigidity [28,29]. Since the bac- 
terial flagellum has a circular cross section, we can choose 
for the material frame the Frenet frame in the undeformed 
ground state, meaning = 0, and thereby restrict the 
spontaneous curvature kq to fl 2 . If and To are con- 
stant along the filament, it has a helical shape with pitch 
p = 2ttto/(kq + Tq) and radius r — ko/( k o + t o)- Note 
that although free energy (3) is formulated in the spirit of 
linear elasticity theory, it becomes highly nonlinear when 
expressed in terms of the space curve r(s) and twist angle 

<t>. 

2.2 Extended KirchhofF rod theory 

In order to describe the transition of the bacterial flagel- 
lum between two polymorphic configurations, the Kirch- 
hoff rod theory has to be extended to include two local 
ground states characterized by fli — (0, Ki,Ti) (i — 1,2). 
We first tried to generalize the approach of Goldstein et 
al., who used a typical double well potential for the twist 
density fl 3 [25,26] (see Appendix A) . However, the force- 
extension curves calculated with this extended free energy 
by numerical simulations did not reproduce the experi- 
mental curves as indicated in Appendix A. 

We, therefore, developed a different model. To each of 
the two relevant polymorphic forms of the flagellum we as- 
sign the elastic free energy (3) of Kirchoff's rod theory and 
introduce a difference S of the two energy densities in the 
ground states. Since the free energy of a system always as- 
sumes a minimum, we locally assign to the bacterial flagel- 
lum with angular strain fl the minimum / po i y (i7, fix, O2) 
of the free energy densities of the two polymorphic con- 
figurations: 

/poi,(«, Ox, fl 2 ) = min (f cl (S2, fl x ), f cl ( fl, fl 2 ) + 5), 

(5) 

where / c i is the elastic free energy density of Eq. (4). The 
resulting density f poly (fl, fix, fl 2 ) is plotted in Fig. 3 as 
a function of curvature k and torsion t using Eq. (2) and 
assuming = 0. 




Fig. 3. (a) Free energy density / po iy from Eq. (5) as a function 
of curvature k and torsion r using fl = (0, k, t). (b) Free en- 
ergy density f polJ as a function of height z and azimuthal angle 
per unit length ip/L of a uniformly stretched helical filament 
(see Section 2.5). The red lines indicate the path of the helical 
filament in the energy landscape during stretching. 

Wada and Netz formulated an alternative, statistical 
model to access the polymorphism of the bacterial flag- 
ellum [27] which they described as a bead-spring chain. 
To each bead they assigned the Kirchoff free energy (4) 
and an Ising spin to distinguish locally between the two 
polymorphic forms of the flagellum. The Ising spin Hamil- 
tonian then favored the same polymorphic state for adja- 
cent beads. By integrating out the spin degree of free- 
dom, they derived an effective free energy density which 
we summarize in Appendix B. In experiments, where a 
thermally induced transition from the normal to the semi- 
coiled configuration was studied, most flagella assumed a 
pure polymorphic form of either the normal or semi-coiled 
state [8] suggesting that the energy cost for forming a do- 
main wall between the two helical states is much larger 
than thermal energy. Using this observation, we demon- 
strate in Appendix B that the free energy density of Wada 
and Netz simplifies to our elastic free energy (5). 

Note that our ansatz can easily be extended to describe 
more than two of the polymorphic states of bacterial flag- 
ella. In particular, this will be necessary for understanding 
the sequence of polymorphic transitions during the tum- 
bling phase of E. coli. 

The elastic free energy density / po i y of Eq. (5) admits 
that two domains of different polymorphic states are sep- 
arated from each other by a sharp domain wall with zero 
width. To realize a more realistic, smoother transition be- 
tween two domains, we have to introduce a free energy 
density of the form [25] 

f bi (d s fi) = ^(d s fi) 2 , (6) 

where 7 1 / 2 is proportional to the width of the domain wall. 
With a properly adjusted 7, experimental observations of 
the domain wall could be described quantitatively [25, 12]. 
Furthermore, our simulations demonstrated that the heli- 
cal filament contains several domains instead of just two 
when 7 is chosen too small or even equal to zero. 

Instead of implementing a constraint for the filament 
to be inextensible, we introduce a stretching free energy 
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density / st = K (d s r) /2. We choose the spring constant 
K such that the changes in the filament length are below 
1.5 %. So the filament is inextensible to a good approxi- 
mation. 

Collecting all the contributions, the total elastic free 
energy T = fds, which we will use in the following for 
modeling the bacterial flagellum, is based on the density 

/ =/poi y («, «i, « 2 ) + f b i(d s O) + f 3t (d s r). (7) 

Since / ultimately depends on the centerline r(s) and 
the twist angle (p(s), the total free energy is a functional 
J 7 [r(s),(j)(s)] in r(s) and 4>(s). 



2.3 Dynamics of a Helical Rod 

We formulate Langevin equations for the location r(s) 
and intrinsic twist <p(s) of the helical filament. At low 
Reynolds number the sum of elastic force per unit length, 
/ el = —SF/Sr, and thermal force / th is balanced by vis- 
cous drag. The same applies to the elastic torque per unit 
length, m e i = —ST jS<\> and thermal torque m t h- Using 
resistive force theory [30] that employs local friction coef- 
ficients 7|i , 7j_ and 7# (see Appendix C), we formulate the 
Langevin equations 

fo|i®i + 7±(l-*®*)]«=/.i + /th (8) 
j r uj =m e i + m th . (9) 

Here t> = dtr is the translational velocity, lo = dt<j> the an- 
gular velocity about the local tangent vector t = e 3 , and 
® means tensorial product. The anisotropic friction tensor 
acting on v in Eq. (8) couples rotation about the helical 
axis to translation and thereby creates the thrust force 
that pushes the bacterium forward [32] . The friction coef- 
ficients per unit length were calculated by Lighthill from 
slender-body theory taking into account the helical geom- 
etry of the rod [31]. The coefficients are summarized in 
Appendix C. In experiments one finds reasonable agree- 
ment with this approach [33,34]. The thermal force / th 
and torque m t h are Gaussian stochastic variables with zero 
mean, (/ th ) = and (m til ) — 0. Their variances obey the 
fluctuation-dissipation theorem and therefore read 

(/«.(*> *) ® / th (*', «')> = 2k B TS(t - t')5(s - a') (10) 
x 7 | 7 1 i®i + 7j 1 (l-i®i) 
(m th (i, s)m th (i', s')) = 2k B TS(t - t')S{s - s')^ 1 



(m th (i, S )/ th (t', S ')> =0. 



2.4 Details of Simulations 



(11) 
(12) 



In our simulations we used a technique developed by Rei- 
chcrt [35] similar to the one of Chirico and Langowski [36] 
and also employed by Wada and Netz in simulations of 
helical nano springs etc. [37,38]. We discretize the center- 
line r(s) of the filament by introducing N + 1 beads at 



locations r % = r(s = i ■ h) and with nearest-neighbor dis- 
tance h. To every bead we attach the material frame, i.e., 
a right-handed tripod of orthonormal vectors {e\, e l 2 , e\} 
(i = 0, . . . , N), where the tangent vector is approximated 

by 



Ti-l 



(13) 



We split the rotation of the tripod along the filament into 
two parts. First, we rotate about the bond direction ef by 
an angle Of h corresponding to the intrinsic twist plus tor- 
sion. Thereafter, the tripod is rotated such that the bond 
orientation e\ is transformed into the consecutive direc- 
tion eg +1 , thus describing the curvature of the filament. 
With this procedure the free energy density / of Eq. (7) 
is discretized and the functional derivatives of the total 
free energy, / el = —ST/Sr and m e i = —ST /Sip, reduce 
to conventional derivatives with respect to r % and (f> 1 . We 
also note that the tangent vector in the friction tensor in 
Eq. (8) is approximated by V = (e| + 4 +1 )/|e| + e* +1 |. 

In the following, we will discuss the influence of the 
three relevant parameters on the force-extension curves; 
the ratio r = C/A of the torsional and bending rigidities 
(twist-to-bend ratio r), the difference in energy S of the 
two helical ground states under study, and the velocity or 
extensional rate v p with which the bacterial flagellum is 
pulled apart. The bending rigidity A together with an ap- 
propriate length introduces characteristic values for force 
and elastic energy. We used A = 3.5pNu-m 2 given in Ref. 
[14] as a typical value for bacterial flagella. All other pa- 
rameters are determined by the geometry of the two poly- 
morphic states. Initially, the bacterial flagellum is in the 
coiled state with spontaneous curvature K\ = 1.8/u.m and 
torsion t\ = 0.56/um and then switches into the normal 
state with k 2 = 1.3/)J.m and t 2 = 2. I/urn. 

The friction coefficients are calculated with the formu- 
las of Lighthill [31] summarized in Appendix C as 7|| = 
1.6-10" 3 pNs/M.m 2 , 7 ± = 2.8 • lO^pNs/nm 2 , and j R = 
0.126 • 10 _3 pNs where a filament diameter of around 20nm 
was used. The length of the filament is L = 10(a.m corre- 
sponding to approximately three helical turns in the coiled 
and four in the normal state. The discretization length be- 
tween the beads was chosen as h = 0.2u-m and the elastic 
coefficient in Eq. (6) as 7 = 0.1pN(j.m 4 . 



2.5 Uniform Deformation 

A helical filament of length L much larger than the pitch 
deforms uniformly aside from inhomogeneities at both ends 
when a constant external force pulls at it. 3 This means 
curvature k and torsion r are constant along the filament. 
We orient the helical axis along the z axis of a cylindrical 
coordinate system (p, ip, z). The geometry and free energy 
of the uniformly deformed helix are completely described 
by its height z — L ^ K i +T i and the difference in azimuthal 



3 For a force compressing the helical filament one observes 
buckling [40]. 



Reinhard Vogel, Holger Stark: Force-extension curves of bacterial flagella 



5 



angles of both ends of the filament, ip = ip(s = L) = 
L^/k 2 + r 2 , where we set ip(s = 0) = 0. Then, the force 
and torque acting on a uniformly stretched helix follows 
from F = —d z T(z,tp) and T = —d v J 7 (z, ip), respectively. 
In experiments, the ends of the filament can freely rotate 
during stretching which means zero torque T. This leads 
to an expression for <p/L 



7_l — (7_l — 7||)(°Wds) 2 . Balancing elastic and viscous 
forces locally, gives the diffusion equation 



ip kq cos a 
L = 



-Tto sin a 



(1 + T) cos 2 a 



(14) 



where we introduced the pitch angle a using the extension 

( = z/L = sin a. (15) 

Under the condition of zero torque, we then obtain the 
classic force-extension relation [28] 

_. . A(rnn cos a + tq sin o)(tq cos a — kq sin a) , „. 

F ^ = — 2 t • 2 , 2 ( 16 ) 

cos a(l cos z a + sin a) z 

_ A ( K o + T o) 3 z a- z ,„ z - z o 
K 2 (r K 2 +r 2 ) L ■ L ' 

In the second line of Eq. (16), the force is linearized 
in a small relative extension (z — zq)/L, where Zq is the 
height of the undcformed helix and k the spring constant. 
Interestingly, we find that curvature k(Q = \J\ — ( 2 <p/L 
and torsion t(£) = (,<p/L of a helix with extension C — z/L 
lie on the ellipse 

[k(0 - «o/2] 2 + r[r(C) - r /2] 2 = (« /2) 2 + r(r /2) 2 , 

(17) 

the center of which is at (no/2, to/2). In deriving Eq. (17) 
we used expression (14) for ip/L. For a symmetric poten- 
tial with twist-to-bcnd ratio r = 1, the ellipse becomes a 
circle. We use relation (17) in Figs. 3(a), 5(b), and 10 to 
indicate the path of a uniformly stretched flagellum in the 
k, t plane. 

Bacterial flagella only have a length of a few pitches. 
So finite size effects from both ends lead to a noticeable 
dependence of curvature and torsion on arc length s. Nev- 
ertheless, we find that Eq. (16) is a good approximation 
for the initial part of the force-extension curve. This is 
valid since the experimental pulling rate is so small that 
the extending filament passes through a sequence of equi- 
librium configurations as we demonstrate in the following 
section. 



2.6 Characteristic Time Scale and Velocity 

Applying instantaneously a force at one end of the heli- 
cal filament induces a localized deformation which spreads 
along the helix. To treat this situation approximately, 
we assume the helical filament to be an extensible rod 
that locally obeys the linearized force-extension relation 
of Eq. (16) with (z — zq)/L replaced by dz/ds. In ad- 
dition, the rod possesses the local friction coefficient per 
unit length for moving a helix parallel to its axis, jh = 



d t z = —d 2 z, 

"fH 



(18) 



where k is the spring constant of the helix given in Eq. 
(16). A localized deformation at one end of the helical fila- 
ment spreads diffusively over the whole filament on a time 
scale re = JhLi 2 /k w 10 _4 s which gives the characteristic 
velocity vc = L/tc ~ 10 3 um/s. In experiments but also 
in our simulations one pulls at the bacterial flagellum with 
velocities much smaller than vc- So the extending filament 
passes through a sequence of equilibrium configurations. 
This also means that the applied and elastic forces nearly 
balance each other since they are much larger than the 
frictional forces. Note that a diffusion equation that bal- 
ances elastic and viscous forces shows up as well in the 
twist dynamics of a rotating elastic filament [44] . 



3 Force-extension curves 

3.1 Pulling on and compressing the helical filament 

We performed both Stokesian (temperature T = 0) and 
Brownian dynamics (T = 300 K ) simulations of the force- 
extension measurements in Ref. [14]. Similar to the ex- 
perimental setup, we fix one end of the filament whereas 
the other end is allowed to move in a harmonic potential 
with trap stiffness 100pN/^m (as in Ref. [14]) mimicking 
the experimental situation where a bead attached to the 
bacterial flagellum was trapped by an optical tweezer. The 
axis of the helical filament is oriented parallel to the z axis. 
In the beginning, the minimum of the harmonic potential 
coincides with one end of the filament in the initial coiled 
state. We then move the potential with a constant veloc- 
ity v p along the z axis and the filament stretches. After 
reaching a maximum extension (m = Zm/L, the potential 
moves with the opposite velocity —v p back into the ini- 
tial position. Note that in experiments the extension rates 
were v p — 0.4u.m/s and less [14]. Such small velocities are 
very time consuming in our simulations so that we typi- 
cally chose values from the range v p = 2, . . . , 20um/s for 
recording the complete force-extension curve. When we 
were just monitoring the initial part of the curve, we used 



velocities as small as v p = 0.2um/s (see Section 3.2). 

We first simulated the extension of a filament with only 
one helical state. Similar to the snapshots 1 and 2 in Fig. 
4(a) the deformation of the filament is uniform except for 
small regions at both ends. This leads to a small deviation 
of the simulated force-extension curve from the theoreti- 
cal prediction of Eq. (16) for a uniformly stretched fila- 
ment. The situation is similar to the results in Fig. 4(b), 
where the initial part 1-2 of the simulated blue curve is 
compared to the analytic prediction (thin black line). If 
we even shift the thin black line to the right, the dashed 
black line agrees very well with the blue curve besides 
the initial part close to position 1. We checked that even 
this difference gradually vanishes with increasing height of 
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(a) (b) 




Fig. 4. (a) Snapshots of a helical filament stretched at T = 
with velocity v p — 2/im/s taken from the supplementary movie 
1. The coiled-to-normal transition (v p > 0) and buckling (v v < 
0) are visible, (b) Force-extension curves simulated without 
thermal noise (T = 0) (blue curve) and with thermal noise: 
two realizations (thin orange and green lines) and an average 
over 10 runs (thick red line) are shown. The parameters are 
v p = 2^im/s and r = 0.7. Analytic prediction for a uniformly 
stretched filament (thin black line) and shifted curve (dashed 
black line). 



the helix as expected. In Section 2.6 we already explained 
that our simulations are performed in the quasistation- 
ary regime. This is also confirmed by our observation that 
at an extension rate v p = 2jj.ni/s the dissipated energy 
during one extension-and-compression cycle is only about 
2% of the maximum elastic energy at extension (m- Fur- 
thermore, we do not observe any pronounced difference 
between deterministic Stokesian and Brownian dynamics 
simulations. For T^0, the forces fluctuate around a value 
which agrees with the deterministic force at T = 0. 

We then determined the force-extension curve at T = 
when the helical filament can switch from the coiled to 
the normal state using the elastic free energy (7). The re- 
sults are illustrated in Fig. 4, where the blue curve in (b) 
corresponds to T = 0. At a certain extension (position 
2) the measured force drops sharply since a small seg- 
ment of the filament close to the fixed end switches into 
the normal state as snapshots 2 and 3 in Fig. 4(a) reveal. 
Then the filament is stretched again. Further adjacent seg- 
ments transform suddenly until at position 6 the filament 
is nearly completely in the normal state. From here we 
invert the velocity v p and move both ends together. How- 
ever, the filament does not transform back into the coiled 
state but remains in the normal state. This ultimately 
leads to a negative force under which the filament buckles 
[see snapshot 7 in Fig. 4(a)]. The full cycle is shown in the 
supplementary movie 1. 

Brownian dynamics simulations reveal the influence of 
thermal fluctuations on the force-extension curve. Two re- 
alizations are included in Fig. 4(b) as thin lines and the 
thick red line shows an ensemble average over 10 runs. A 
full cycle of one realization is shown in the supplemen- 



tary movie 2. Clearly, the first transition into the normal 
state occurs at a smaller extension compared to the deter- 
ministic case since thermal fluctuations help to overcome 
the energy barrier in the elastic free energy density (see 
Fig. 3). Whereas the force in each realization fluctuates 
visibly, the sharp decrease of the force when a local tran- 
sition into the normal state occurs is as pronounced as in 
the deterministic case. Finally, when both ends of the fil- 
ament are moved together, it completely transforms back 
to the coiled state and buckling is not observed. The sin- 
gle realizations of the complete cycle of the force-extension 
curve closely resemble the experimental curves in Figs. 5 
and 6 of Ref. [14]. In particular, the force and extension 
where the first coiled-to-normal transition occurs fall into 
the experimental ranges of 3 to 5 pN and ( = 0.55 to 0.6, 
respectively. 

We now discuss in more detail how the difference 6 
in the ground-state energies of the two helical states, the 
twist-bend ratio F, and the extension rate v p influence the 
force-extension curves. 



3.1.1 Ground-state energy difference 5 of the coiled and 
normal state 

Increasing the ground-state energy 6 > of the normal 
relative to the coiled state also increases the energy bar- 
rier, which the flagellum in the coiled configuration has 
to overcome to transform locally into the normal state. 
Therefore, the transition is delayed to a larger extension 
z/L or does not occur at all. On the other hand, the bar- 
rier which the normal configuration has to overcome to 
relax back into the coiled state decreases and buckling 
of the filament becomes less probable. Observations also 
show that a filament prepared in the normal state very 
slowly relaxes back into the coiled state, so 5 should not 
be too large compared to thermal energy ksT. Using also 
the following quantitative considerations, we adjusted S to 
O.lpN which resulted in the good agreement with experi- 
mental observations, already demonstrated. 

We now derive an upper bound for S to ensure that a 
transformation from the coiled to the normal state is, in 
principal, observable. The locus of the energy barrier in 
the K, t plane of Fig. 3 is given by 

f cl (n,n 1 ) = f cl (n,n 2 ) + s, (19) 

where Q = (0, k,t) and = (0, Kj,Tj) contains the val- 
ues for spontaneous bend Ki and torsion Tj of the coiled 
(i = 1) and normal state (i = 2), respectively. Assum- 
ing both states have the same elastic constants A, C, we 
arrive at a straight line 

k(«i - n 2 ) + tT(ti - t 2 ) + 5/ A 

= ( K 2 1 -4)/2 + r(T 2 1 -T 2 2 )/2. (20) 

We will use this formula later. The ground-state energy 
density of the normal state at fi 2 = (0, «2,Ta) is S. Only 
if S lies below the energy density of the coiled state at 
fl = fl 2 is a clear transition between both states possible. 
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The explicit form of this upper bound, 6 < f c i(f2 2 , 
combines with S > to the inequality 



< 



A(K! -k 2 ) 2 



<(i + rz\ 2 )/2, 



where 



n - r 2 

Ki - K 2 



(21) 



(22) 



is the ratio of the differences in spontaneous torsion and 
curvature. 

In experiments the value of 5 changes with the condi- 
tions of the solvent. In particular, different polymorphic 
forms of the bacterial flagcllum become stable when one 
alters the pH value, ionic strength, or temperature of the 
aqueous environment [8]. It would be interesting to per- 
form the force-extension experiments under different con- 
ditions to investigate how a changing 5 but also variations 
in bending (A) and torsional (C) rigidity influence the 
force-extension curve of a bacterial fiagellum. 



3.1.2 Twist-to-bend ratio T 

Increasing the twist-to-bend ratio r also increases the en- 
ergy barrier as the following formula for the minimum 
value of the barrier demonstrates, 



fb Jl + 25/[A( Kl -K 2 ) 2 ]+rA 2 
A{k 1 -k 2 ) 2 8(1 + /Y\ 2 ) 



(23) 



However, in contrast to S an increasing _T increases both 
barriers for the coiled-to-normal and the normal-to-coiled 
transition. We now discuss in Fig. 5(a) how the twist-to- 
bend ratio r influences the force-extension curve of the 
helical filament. The filament is stretched to an extension 
( = 0.95, well above the equilibrium height ( = 0.85 of 
the normal state. Fig. 5(b) shows contour plots of the 
elastic free energy density as a function of curvature and 
torsion for the same r as in (a). The red line indicates 
the sequence of k, t values, as the filament in the coiled 
state is uniformly stretched (see Section 2.5). Note that 
for r ^ 1 the anisotropy of the elastic free energy density 
is clearly visible. 

For r — 0.5 and 0.7, one observes the typical force- 
extension curves as already discussed in Fig. 4 also for 
r = 0.7. Since now the maximal extension ( = 0.95 is 
larger, the curve for r = 0.7 during compression looks 
different. While at £ = 0.95 the entire filament is in the 
normal state, the filament for r = 1 no longer transforms 
completely into the normal state. Most pronounced at 
r = 1 is the fact that even with thermal forces the filament 
does not return into the coiled state during compression. 
This results in a negative force under which the filament 
buckles. All the force-extension curves in Fig. 5(a) are de- 
termined for an extension rate v p = 2um/s. A smaller v p 
enhances the probability that thermal fluctuations trans- 
form the filament back into the coiled state and buckling 
is not observed. The same is true for a smaller extension 
where a larger part of the filament stays in the coiled state 



and makes it easier for the rest of the filament to return 
to the coiled state. At a ratio r — 1.5 the force-extension 
curve changes qualitatively. A first transition into the nor- 
mal state occurs at a larger extension £ w 0.65 compared 
to the previous curves. The transitions are no longer as 
pronounced and there are larger differences between sin- 
gle realizations of the force-extension curve. In addition, 
for r < 1.5 only one normal domain occurs whereas for 
r = 1.5 two normal domains are observed. The reason for 
all these features becomes clear from Fig. 5(b). At r = 1.5 
a uniformly stretched filament does not hit the energy bar- 
rier anymore but passes the barrier in a close distance so 
that thermal fluctuations are needed to induce transfor- 
mations into the normal state. For further increase of JT 
a transition into the normal state does not occur at all 
realizations, as the three graphs for r — 2 demonstrate. 
A similar behavior is observed in Ref. [14] for a transition 
from the normal into the hyperextended state. The blue 
curve corresponds to the traditional force-extension rela- 
tion of the helix in the coiled state. In the yellow curve 
a complete transition into the normal state was realized. 
A partial transformation occurred in the magenta curve 
which then completely returned into the coiled state. Fi- 
nally, at r = 2.5 the filament always remains in the coiled 
state. Beyond an extension of £ = 0.5 the slope of the 
curve becomes smaller. This is the onset of a qualitatively 
new behavior. At even larger r and t /k < 1, one ob- 
serves a sharp drop in the force due to a discontinuous 
transformation, where one turn of the helical filament un- 
winds. This is discussed in Refs. [39,37]. 

As discussed, neglecting thermal fluctuations, a coiled- 
to-normal transition is no longer observable in the force- 
extension curve when in Fig. 5(b) the straight line [Eq. 
(20)] separating the coiled and normal state becomes tan- 
gential to the trajectory [Eq. (17)] of the uniformly stretched 
filament. Combining both equations leads to a second con- 
dition for S I A: 



< 



1 



A(ki - k 2 ) 2 2(«i - K 2 ) 



K2 + rr 2 A 



^l + rA 2 )(K 2 +rr 2 )) (24) 



Together with condition (21), we obtain a region in the pa- 
rameter space (r,5/A) where a transformation from the 
coiled to the normal state should occur. The region is in- 
dicated as shaded area in Fig. 6. Based on experimental 
data for Young's and shear modulus in literature, Flynn 
and Ma received for the twist-to-bend ratio r the range 
2 • 10 _1 -2 • 10 2 [41]. With a computational method, called 
the quantized clastic deformational model, they calculated 
r ss 23 [41]. On the other hand, together with our value 
5/[A{m — k 2 ) 2 } « 0.11, we predict a value of r < 1, in 
good agreement with our simulations. 



3.1.3 Extensional rate 

In Section 2.6 we reasoned that during the measurements 
of the force-extension curves the bacterial fiagellum goes 
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Fig. 5. (a) One or several realizations of force-extension curves for increasing twist-to-bend ratios r. The extensional rate is 
v p = 2/im/s. (b) Contour plots of the elastic free energy density as a function of curvature and torsion for the same F as in (a). 
The red line indicates the sequence of k, t values, as the filament in the coiled state is uniformly stretched. 





Fig. 6. The shaded area indicates the parameter ranges for the (c ^ d) U . J . u '' (n u nmal) (foiled) ^ ■ °' 7 (normal 

ground-state energy difference 8 and the twist-to-bend ratio extension £ extension £ 
r for which a coiled-to-normal transition should be observed. 

The dots give parameter values used in simulations. Note that Fig. 7. Force-extension curves for twist-to-bend ratios r = 0.7 

thermal fluctuations shift the upper border for F indicated by (a) and r = 1 (b) at different extension rates v p = 2u.m/s, 



the blue line to the right. 



5u.ni/s and 20um/s. 



through a sequence of equilibrium states. This means that 
frictional forces acting on the filament through the sol- 
vent are negligible against elastic forces. Therefore, with- 
out thermal noise, the force-extension curve does not de- 
pend on the extensional rate v p . In contrast, our Brown- 
ian dynamics simulations demonstrate a clear influence of 



v p . In Fig. 7, we show force-extension curves for twist-to- 
bend ratios r = 0.7 and 1 and different extension rates 
v p = 2u.m/s, 5u.m/s and 20u.m/s. The first transition from 
the coiled to normal state occurs at larger extensions when 
v p is increased. This is immediately clear since smaller ve- 
locities v p give the filament more time to explore the en- 
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Fig. 8. Mean-first-passage time (MFPT) as a function of ex- 
tension £ for f = 0.7 and _T = 1.0. Inset: When the local 
coiled-to-normal transition occurs, the energy decreases instan- 
taneously. 
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ergy landscape with the help of thermal fluctuations. So 
the appropriate local curvature and torsion to overcome 
the potential barrier can be created at smaller extension. 
The curves in Fig. 7 just give one specific realization for 
each parameter set. In the following section, we will inves- 
tigate in detail the probability distribution for the exten- 
sion where the first coiled-to-normal transition occurs. 

If the filament in the normal state is compressed too 
fast, it will start to buckle since again it has not sufficient 
time to overcome the energy barrier. This is visible in 
Fig. 7(a), where the filament with the highest compressing 
rate v p — 20/xm/s buckles which corresponds to a negative 
force. On the other hand, for smaller rates v p , the filament 
always returns into the coiled state. 



3.2 Clamped filament 

Since the bacterial flagellum goes through a sequence of 
equilibrium states during the measurements of the force- 
extension curves, it should be possible to derive some char- 
acteristics of these curves by investigating a clamped fila- 
ment which is hold at a fixed extension £ = z/L in thermal 
equilibrium. In particular, we show here how one can in- 
fer the mean extension (Q(v p ) at which the filament in 
a force-extension measurement undergoes the first local 
transition to the normal state 

We stretch the filament in the coiled state at zero tem- 
perature to an extension £ and keep this extension con- 
stant. We then perform a Brownian dynamics simulation 
and determine for each realization the first-passage time 
r at which a local transition to the normal state occurs. 
The inset of Fig. 8 shows the filament before and after the 
transition accompanied by a sharp decrease of the elastic 
free energy. According to Kramers theory, the mean-first- 
passage time (MFPT) is proportional to the Arrhenius 
factor, 



ib" 1 10° 10 1 10 2 

extensional rate v p [^m/s] 

Fig. 9. (a) Probability distributions p(£, v p ) as a function of 
extension £ for r — 0.7 (blue) and JT = 1 (magenta) at exten- 
sion rates v p — 0.2fi.m/s (full lines) and v v = 20u.m/s (dashed 
lines), (b) Analytically determined mean extension {£) (v p ) for 
the coiled-to-normal transition as a function of velocity v p for 
r = 0.7 (blue) and r = 1 (magenta). The dots indicate nu- 
merically determined extensions for each simulation run. 

The activation energy AT(() depends on the extension £ 
and is needed to create the domain wall between the coiled 
and normal state. We determined the MFPT by averaging 
over 100 simulated values for r at each extension £. The 
results are plotted in Fig. 8 for two twist-to-bend ratios 
r = 0.7 and r = 1. We only simulated (r) for a small 
variation in the extension for two reasons: at larger ex- 
tensions where A.F(£) ps fc B T the Kramers theory is no 
longer valid and (r) is too small to be determined accu- 
rately; at smaller extensions (t) is so large that it cannot 
be calculated in reasonable simulation times. Note that 
for r = 1 the MFPT spans three decades. Our simulation 
results demonstrate that log (r) (() can be fitted by 

log (r) (0 w a - PC, (26) 

where (r) is measured in seconds. Equation (26) can just 
be viewed as a Taylor expansion to linear order in £. The 
parameters a and f3 follow from a least-square fit. Whereas 
a changes with r, the slope (3 surprisingly does not seem 
to depend on _T within the numerical accuracy. A similar 
law as Eq. (26) is used in situations where bonds rupture 
under a given load force [43,45]. Note, however, that here 
we control the extension. 

Within the adiabatic or quasistationary approxima- 
tion, we now formulate the probability p(t, v p )dt that with- 
in the time interval [t, t + dt] the filament transforms lo- 
cally from the coiled to normal state when stretched with 
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velocity v p , 

Pit, v p ) = ex P (- f Q 1/ <r> («t',v P W) , 

(27) 

where 

{{t,v p ) = Co + tv p /L (28) 

is the extension at time t. The exponential factor in Eq. 
(27) is the probability that until time t a transition does 
not occur, which follows by writing the probability that 
a transition within the time interval dt' does not hap- 
pen as (1 — dt' / (t)) ?s exp(— dt' / (r)). The first factor is 
the probability per unit time that the transition occurs 
within dt at time t when the filament is with certainty 
in the coiled state before. Using Eq. (28), one introduces 
the probability p((,v p )d( = p(t,v p )dt that the filament 
undergoes a local coiled-to-normal transition at extension 
(. We calculate it analytically in Appendix D by assum- 
ing the validity of Eq. (26) for the whole £ range. The 
results are plotted in Fig. 9(a) for r = 0.7 (blue) and 
r = 1 (magenta) and for extension rates v p = 0.2pm/s 
(full lines) and v p = 20[J.m/s (dashed lines). The probabil- 
ity distributions are concentrated on a small range about 
their maximum values and are shifted to larger extensions 
for increasing velocities v p , as expected. We then calculate 
the mean extension 

/>oo 

(C)K)= / Cp(C,v P )dC, (29) 

"'Co 

at which the coiled-normal transition occurs first for a 
given extension rate v p . The evaluation of the integral is 
performed in Appendix D and finally gives 

(C) (v P ) oc logv p . (30) 

The mean extension (£) (v p ) with all prefactors and con- 
stant terms is plotted in Fig. 9(b) for r = 0.7 (blue) and 
r = 1 (magenta) together with numerically determined 
extensions £ for several realizations at a specific velocity 
v p . The values scatter around the mean value and are, 
therefore, in good agreement with the analytical treat- 
ment. We note that a relation similar to Eq. (30) occurs 
for the velocity dependence of the rupture force of single 
molecular bonds in dynamic force spectroscopy [43]. 

4 Summary and conclusions 

In this article we have developed a sufficiently simple elas- 
tic model to describe the polymorphism of a bacterial 
flagellum based on Kirchhoff's theory of an elastic rod. 
The friction with the aqueous environment is modeled 
within resistive force theory. Using geometrical parame- 
ters of the coiled and normal states and the bending rigid- 
ity as obtained in Ref. [14], we are able to reproduce the 
force-extension curves recorded in experiments. Thermal 
fluctuations realized within Brownian dynamics simula- 
tions are crucial. In particular, the force values at which 



a first coiled-to-normal transition takes place lie between 
3 and 5 pN, as in experiments [14]. We have investigated 
in detail how the force-extension curve depends on the 
twist-to-bend ratio r and, furthermore, by analytic argu- 
ments identified a parameter region for ground-state en- 
ergy difference S and -T, where a coiled-to-normal transi- 
tion should be observable. It clearly demonstrates that for 
values of r well above one, a polymorphic transformation 
is not possible and therefore contradicts some of the ex- 
perimental values for r recorded in literature. Based on 
our simulations, we predict r to be within 0.7 and 1.0. 
Further studies demonstrate how the extensional rate v p 
influences the force extension curve. We directly observe 
the influence in the simulations of the full force-extension 
cycle but also when we concentrate on the extension for 
the first coiled-to-normal transition. Since the extensional 
rates v p are sufficiently small, the flagellum goes through 
a sequence of quasi-stationary states. We, therefore, used 
equilibrium properties of clamped flagella to predict a 
logarithmic velocity dependence for the mean extension 
of the first coiled-to-normal transition in good agreement 
with our simulations. 

Our approach is easily extended to include more than 
two polymorphic states. Figure 10 shows the contour lines 
of the resulting free energy landscape when all twelve poly- 
morphic conformations of Fig. 1 are taken into account. 
For each of these conformations we take Kirchhoff's elas- 
tic free energy with spontaneous curvature and torsion as 
given by Calladine [17] and just consider the minimum 
value from all these free energies. We assume here that 
the ground-state energies of all helical states are zero and 
that they all have the same bending and torsional rigidi- 
ties with r = 0.7. We also included the paths of uniformly 
stretched helical filaments. The coiled-to-normal transi- 
tion (nn — 3 to 2) takes place with certainty. However, 
all the other transitions would need thermal fluctuations 
to occur. For example, the normal-to-hyperextended tran- 
sition (riR = 2 to 1) will only occur when it is stretched 
sufficiently slowly so that thermal fluctuations can induce 
the transformation. 

We can now use our model to study various aspects 
connected to bacterial locomotion. For example, we will 
investigate how polymorphic transitions are induced by 
rotating the flagellum. Our very challenging goal is to 
model the complete tumbling cycle of a bacterium. Fi- 
nally, we note that our model may be applicable as well 
to spirochetes, where very recent experiments and theory 
have begun to address the elasticity of the helical structure 
[46]. 
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Fig. 10. Contour lines of the free energy density of the extended Kirchhoff rod theory that includes all 12 polymorphic states 
of the bacterial flagelhim. For all helical states the same ground state energy and the same bending and torsional rigidity with 
r — 0.7 is assumed. The red lines indicate the paths of uniformly stretched flagella. 



A Extended Kirchhoff rod theory: the 
double-well potential 

Here, the main idea is to extend the theory of Goldstein 
et al. [25,26]. They used a double well potential for the 
twist density Q3, realized by a polynomial of degree four, 
to describe two polymorphic states of a flagellum [25,26]. 

In order to develop a strategy how to generalize this 
double well potential to all three coordinates Qt (i — 
1, 2, 3), we first write down a general one-dimensional poly- 
nomial of degree four: 



f{x,x l7 x 2 ) 



A (x-xif 



2 {x 1 
_ d 
~ 6 



x 2 y 



\(x-x 2 f 



(x-ari)(3a: + xi-4!Ca)]. (31) 



For d < 1 it has two minima at x\ and x 2 with f{x\) 
and f{x%) — 5, respectively, where 







S = J^ A (xi - x 2 fd. 



Whereas d^f(xi) — A, the second derivative at x 2 de- 
pends on the parameters x±, x 2 , d, and A. 

We generalize the polynomial of Eq. (31) to three di- 
mensions by replacing terms of the form Axy by x ■ Ay, 
where a?, y are three-dimensional vectors and A is a di- 
agonal matrix with An — A 22 — A and A33 = C. The 
constants A, C are the bending and torsional rigidities, 
respectively. Using the shorthand notation | 
with x = {2 — fli. we now write 



f(f2,f2 u f2 2 ) = 



1 \n 



2\f2 1 -f2 2 \i' 



\n- 



-(O - Oi) • A(3f2 + f2i- 4J? 2 )] 



1 

12' 



with S = -^1^1 - n 2 \\d. 



(33) 
(34) 



The polynomial is illustrated in Fig. 11(a) with fl = 
(0,k,t). It has two minima at fl\, fl 2 with f{fl\) = 0, 




0.3 
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0.7 0.85 
extension ( (normal) 



Fig. 11. (a) Double- well potential of Eq. (33) plotted as a 
function of Q2 = k, and O3 — r for the parameters of the 
coiled and normal helical state, (b) Force-extension curve (red 
line) simulated with this potential at T — and v p = 2(i.m/s. 
The blue curve is taken from Fig. 5(a) for the same parameter 

r = 1. 



f{fl 2 ) — S, respectively, and one saddle point at f2 3 
(32) fii + 2^7;(^2 — ^1) for d < 1. Close to the first minimum 



at J?i, the polynomial agrees with Kirchhoff's elastic free 
energy (4) , whereas the bending and torsional rigidities at 
fl 2 and also the energy barrier depend on 8. 

Figure 11(b) shows a force-extension curve (red line) 
simulated with Eq. (33) at T = for the same parame- 
ters as in Fig. 5(a) with r = 1. For comparison the initial 
part of the curve from Fig. 5(a) (r = 1) is included. The 
helical filament is much softer and the initial part of the 
force-extension relation has a negative curvature in con- 
trast to experiments. We, therefore, decided to introduce 
the alternative model of Eq. (5). 



B Simplifying the free energy of Wada and 
Netz 

After integrating out the spin degree of freedom, Wada 
and Netz arrived at the following elastic free energy den- 
sity [27]: 



(35) 
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with 

h = 



K\ + K 2 



c 



(36) 



r kBT ^ 
I 2 = In 



cosh(yl) + Jsinh 2 (A) + e - iJ / k ^ T 



J 

a 



(37) 



and 
k B T 



a r & \ ( n Ki+ K 2 

-A=S + -{K 1 -Ka)(fi2-— r - 



+ ^(ri-r 2 ) ( il. 



(38) 



Here A, C are the bending and torsional rigidities, respec- 
tively, 5 is the difference of the ground-state energies of 
the helical conformations (ki,ti) and (k 2 ,t 2 ), and a the 
length of discretization. The quantity J is the interaction 
strength in the Ising Hamiltonian. 

We interpret the energy cost 23 for two anti-parallel 
spins as the energy of a domain wall connecting two he- 
lical states. Since such domain walls are rarely seen in 
experiments [8], we can assume 2J >• k B T and, therefore, 
approximate f 2 as 



k -^\A\ + J -, 
a a 



(39) 



where we have used cosh a; + | sinhx| = exp \ x\. Now, we 
introduce the elastic free energy densities of the two helical 
states, 

a = + ^ {Q2 ~ K2)2 + ^ {f2s ~ T2)2 + S/2 > (40) 

/3 = \n\ + ^{n 2 - Kl ) 2 + ^ (n 3 - Tl ) 2 - 6/2, (4i) 

and write the two contributions to the energy density (35) 
as 



/2«|l«-/?l + ^» (42) 



where c\ is just a constant. This finally gives our ansatz 
for the free energy density (5): 

f^(a + p-\a-p\) + Cl =min(a,/3) + Cl . (43) 



C Friction coefficients per unit length 

In a moving helical filament, different parts interact via 
hydrodynamic interactions. Nevertheless, using slender- 
body theory, Lighthill demonstrated that one can describe 
the hydrodynamic friction of the filament with the help of 
resistive force theory that introduces local friction coeffi- 
cients per unit length. For the translational coefficients, 
Lighthill obtained [31] 



7|| 



2tti] 



\n{2q/r) 



and 7j_ 



ln(2g/r) + 1/2 



(44) 



Here rj is the shear viscosity, r = 0.02ixm the cross-sectional 
radius of the bacterial flagellum, and q a characteristic 
lengt h, for wh ich Lighthill derived q = 0.09^1, where A = 
27t/v / k 2 + t 2 is the filament length of one helical turn. 
This gives 7j| « 1.54/7, 7j_ « 2.747? for the normal state 
and 7|| ~ 1.647;, 7j_ w 2.91r) for the coiled state. Since the 
coefficients arc similar in both states, we use the interme- 
diate values 7|| ~ 1.677 an d 7-L ~ 2.877 f° r simulating the 
force-extension curves. Finally, for the rotational friction 
coefficient one finds [42] 

7 fl = 47r/ja 2 . (45) 

D Velocity dependence of the mean 
extension (() (v p ) 

In the following we use the numerically determined MFPT 
of Eq. (26) in the form 



(r) =r cxp[-/3(C-Co)]. 



(46) 



We transform the probability distribution p(t) in Eq. (27) 
into a probability distribution for the rescaled extension 
variable x — (3(( — (0) = (3v p t/L using p(t)dt = p(x)dx. 
The integral in Eq. (27) is readily calculated with the help 
of (t) = r exp(— x) and we obtain 



p(x)dx = e x A exp (-e A (e x - 1)) da:, 



(47) 



where exp(— A) = L/((3v p t ). At x = or ( = Co, the 
choice of our parameters shows that the energy barrier 
between coiled and normal state is much larger than k B T, 
so p(x = 0) = e~ A is almost zero or A> 1. Even for x 
0, we can therefore approximate the probability density in 
Eq. (47) as 



p(x) w c x ' A exp (-e x A ) Ax := p a (x - A) 



(48) 



It explains why all the probability densities p(C,,vp) in 
Fig. 9 have the same shape but are shifted relative to 
each other along the £ axis due to different values of v p 
and r and, therefore, of A. The distribution p(x) has a 
maximum at x = A with p(A) we -1 . 

We now calculate the rescaled mean extension 



(x 



/>oo 

) = / xp(x)dx. 
Jo 



(49) 



We rewrite the probability distribution of Eq. (48) as 
p(x) = c x ~ A exp (— c x ~ A ) = —d x exp (— c x ~ A ) and obtain 



after integrating by parts 



(x) = cxp(e ) / exp (- 
Jo 



)dx. 



(50) 



The substitution y — — exp(x — A) then introduces the 
exponential integral function: 



(x) — — exp(e 



(51) 
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which for e A <C 1 is approximated by 

(x) «-(C + logc- A ) = -C + A, (52) 

where C s=s 0.577 is Euler's constant. Introducing the orig- 
inal extension variable ( and A = \og(v p T n /L) + log(/3), 
we obtain for the mean extension at which the coiled-to- 
normal transition first takes place: 

/3«C) - Co) = -C + log(v p T Q /L) + log(/?) cx log(v p ). 

(53) 
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